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Abstract. We derive improved versions of the relations between supermassive black hole mass 
(Mbh) and host-galaxy bulge velocity dispersion (<r) and luminosity (L) (the M-a and M-L 
relations), based on ~ 50 Mbh measurements and ~ 20 upper limits. Particular attention is paid 
to recovery of the intrinsic scatter (to) in both relations. We find the scatter to be significantly 
larger than estimated in most previous studies. The large scatter requires revision of the local 
black hole mass function, and it implies that there may be substantial selection bias in studies 
of the evolution of the M-a and M-L relations. When only considering ellipticals, the scatter 
decreases. These results appear to be insensitive to a wide range of assumptions about the 
measurement errors and the distribution of intrinsic scatter. We also investigate the effects on 
the fits of culling the sample according to the resolution of the black hole's sphere of influence. 
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1. Overview 

The M-L and M-a relations — the relations between a black hole's mass and the host 
galaxy's (bulge) luminosity, L, or velocity dispersion, a — strongly suggest a fundamen- 
tal link between galaxy and black hole (BH) evolution (iDresslerl 19891 : Kormendv 1993 : 
iMaeorrian et al.lll998HGebhardt et al.ll2000tlFerrarese fe Merrittll20obT ). In this contribu- 
tion to the proceedings, we discuss recent developments in the study of these relations 
and how they relate to coevolution of black holes and galaxies. Fundamental to the un- 
derstanding of the M—a and M-L relations is the measurement of the relation's intrinsic 
or cosmic scatter, as distinct from scatter due to measurement errors. The fact that there 
is a relation between BH mass and stellar velocity dispers i on is n ot surprising, but the 
scatter is remarkably small, estimated bv iTremaine et all ( 2002 ) to be no larger than 
0.25-0.3 dex. The total scatter of the relations (intrinsic combined with statistical and 
systematic measurement errors) is what makes the relation useful as a secondary tool for 
BH mass estimation. The intrinsic scatter, however, is the measure of the fundamental 
link between the physical quantities in question. 

The magnitude of the intrinsic scatter is extremely important for several reasons. First, 
the range of BH masses in galaxies of a given velocity dispersion or bulge luminosity 
constrains BH formation and evolution theories. Many theories of BH formation and 
galaxy evolution have used the M-a relation either as a starting point for furthe r work 
or as a pr ediction of the th eory (e.g., Silk fc Rees 1998 : Fabian 1999t lKmg]|2003 ); for a 



review, see lRichstond(|2004l ). A further test of such theories is whether they can reproduce 



the observed cosmic scatter in the relation. For example, there may be an increased 
intrinsic scatter in low-mass galaxies because BHs are ejected by asymm etric gravi t ationa l 
wave emission and lo w- mass spheroids have lower escape velocities IVolonteril (|2007l ); 
Volonteri et "all (|2008f) . 
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Understanding the scatter in the M-a relation is also essential for estimating the space 
density of the most massive BHs in the local universe. One of the most useful aspects 
of the M-a relation is that it allows one to estimate a galaxy's central BH mass from 
the more easily measured velocity dispersion. Becau se of the steep de c line in number 



density of galaxies hav ing high velocity dispersion (jSheth et al. 1 120031: iBernardi et al 



20061: lLauer et al. I l2007h . the majority of the extremely large BHs will reside in galaxies 
with moderate velocity dispersi ons that happen to contain BHs that are overmassive 
for th e given velocity dispersion (|Yu fc TremaindEooi Im arconi et al. 2004: Lauer et al 



20071 ) . Knowing the magnitude of the intrinsic scatter is thus required to find the density 
of the most massive BHs. For example, the number density of BHs with M > 10 10 M e 
3 Gpc ~ 3 if the intrinsic scatter is 0.15 dex and ~ 30 Gpc~ 3 if the intrinsic scatter 



is 



is 0.30 dex ( Lauer et al.ll2007h 



Both the magnitude of the intrinsic scatter and its distribution (e.g., normal or log- 
normal in ma ss) are also important to know for studies of the evolution of the M-a 
relation (e . g.. iTreu et al.l I2004L 120071: iHopkins et "all 120061 : IPeng et all 120061 : IShen et al ' 



2007tl2008tivestergaard et alJl2008h . lLauer et all (|2007h ~ showed that there is a bias when 



comparing BH masses derived from observations of inactive galaxies at low redshifts to 
BH masses from active galaxies at higher redshift. The bias arises because the sample 
of nearby galaxies measures the distribution of BH masses for a given host velocity 
dispersion or luminosity, whereas the sample from high-redshift galaxies tends to measure 
the distribu t ion of the host luminosity or host velocity dispersion for a given BH mass. 
Lauer et al. I (l2007h found that the bias in the inferred logarithmic mass scales as the 
square of the intrinsic scatter in logarithmic mass. In order to account for this bias 
correctly, not only the magnitude but also the distribution of the deviations from the 
M-a relation is needed. 



2. Most Recent Scaling Relations and their Scatter 



Using a sample of ~ 50 BH mass measurements and ~ 2Q. lGiiltekin et al.l (|2009bT ) used 
a generalized maximum likelihood method to fit the M-a and M-L relations with an 
intrinsic scatter component. Figure [T] and figure [2] show M-a and M-L relations. Their 
best fit for M-a was 

log (^) = 812 ± °- 08 + (4 ' 24 ± a41) log (2601^=1) (2 ' 1} 

with an intrinsic scatter distributed normally in logarithmic mass with standard deviation 
eo = 0.44 ± 0.06. When considering only ellipticals, they found that the intrinsic scatter 
drops to eo = 0.31 ± 0.06. Their best fit for M-a was 



M BH \ 



log [j^j = 8-95 ± 0.11 + (1.11 ± 0.18) log [j^J (2-2) 

with an intrinsic scatter distributed normally in logarithmic mass with standard deviation 
eo = 0.38 ± 0.09. Different assumptions about the error distribution in black hole mass 
measurements and intrinsic scatter did not significantly alter these conclusions. 



3. Implications for Scatter: Potential for Bias in Co-evolution Studies 

The intrinsic scatter is essential for determining the cosmic density of the most massive 
black holes. Because of the exponential drop in number density at the high end, very 
large galaxies are extremely rare. This means that in a fixed volume, the greatest number 




Figure 1. The M—a relation for galaxies with dynamical measurements. The symbol indicates 
the method of BH mass measurement: stellar dynamical {pentagrams), gas dynamical {circles), 
masers {asterisks). Arrows indicate 3<76g upper limits to BH mass. The shade of the error ellipse 
indicates the Hubble type of the host. The saturation of the shades in the error ellipses is 
inversely proportional to the area of the ellipse. The line is the best fit relation to the full 
sample: Af B H = 10 812 M o (cr/200 km s" 1 ) 4 24 . The mass uncertainty for NGC 4258 has been 
plotted much larger than its actual value so that it will show on this plot. (Adapted from 
iGiiltekin et~alll2009bl .') 



of big black holes do not come from the intrinsically-rare, large galaxies but from the 
more common modest-sized galaxies that happen to have an over-massive black hole. 
Thus, deriving the number density of the largest black holes from the number density of 
galaxies depends on this effect. 

The tendency for the largest black holes to come from more modestly si zed galaxies 



also l eads to a potential bias in studies of the evolution of scaling relations (jLauer et al 



20071) . The samples of high-redshift black holes, which are measured from AGNs, tend 



to probe distribution of a or L for a given black hole mass whereas the local, quiescent 
sample tends to probe the distribution of Mbh for a given host galaxy property. Direct 
comparison of the two samples leads to a bias that would lead one to incorrectly infer 
that at high rcdshift black holes were more massive than they are at low redshift, for a 
given velocity dispersion or bulge luminosity. This bias depends critically on the poorly 
constrained wings of the distribution of the intrinsic scatter. 

Another, previously neglected caveat is that it may be impossible to disentangle evo- 
lution of just the intrinsic scatter of the relations from the evolution of the slope and/or 
intercept of the relations. That is, to correct for this effect when looking for evolution 



122 



Kayhan Giiltekin 



10" 



n — i — i — | — I — i — I — i — | — I — i — r 

A Stellar Kinematics 
O Gas KiaeawlKS 

• Elliptical 

• SO 



-i — i — r — l — I — r — l 




IB* 



'Ml 



0-5 



lfljO 



L1.0 



Figure 2. The M-L relation for galaxies with dynamical measurements. The symbol indicates 
the method of BH mass measurement: stellar dynamical (pentagrams) and gas dynamical (cir- 
cles). Arrows indicate upper limits for BH mass. The shade of the error ellipse indicates the 
Hubble type of the host galaxy, and the saturation of the shade is inversely proportional to 
the area of the ellipse. The line is the best-fit relati on for the sample wit hout upper limits: 
M BH = 10 8 95 M o (L v /10 n Lq.v) 111 . (Adapted from lGiiltekin et alj|2009bh 



of M-a or M-L, one needs to assume that at high redshift the intrinsic scatter is the 
same as at low redshift in the quiescent population. This assumption that one aspect 
of the relation (the scatter) is constant while measuring changes in other aspects (the 
intercept and slope) i s difficult to justify a nd may lead to erroneous conclusions. Finally, 
a very recent study ( Shen fe Kelly 20091 ) has shown that an independent bias arising 
from uncertainties in BH mass estimators and the shape of the BH mass function may 
lead to an overestimate of the true BH masses. 



4. Scaling Relations are Biased When Black Hole Masses are 
Censored by Sphere-of-Influence Resolution 

Some have advocated the removal of BHs from scaling relati on estimation if the BH 's 
radius of influence on the sky is below a given threshold (e.g., iFerrarese fe Ford 2005 ). 
This is based on a misconception that under-resolved BHs will yield biased BH mass 
estimates. In fact, under-resolve d BHs yield lower precision masses, but with no bias in 



accuracy (|Giiltekin et al.ll2009alfbh . On the other hand, we may demonstrate the effects 



of censoring data based on i?; n fiwith a synthetic M-a data set. The synthetic data set 



Scatter in M-a and M-L 
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Figure 3. This figure shows results of a synthetically generated sample of 500 galaxies with 
BH mass generated from an M-a relation with a — 8 and j3 = 4.0 and a log-normal scatter 
of 0.3 dex with measurement errors of 0.2 dex. The M-a ridge line is drawn as a black line. 
The Galaxy is plotted as a pentagram. Different symbols indicate different levels of resolution: 
ftnfi/j" rC s < 1-0 (squares), 1.0 < i?i fl/r r8 s < 2.0 (diamonds), 2.0 < Ri n a/r TC s < 4.0 (triangles), 
and -Rinfl/r'rcs > 4.0 (circles). Fitting the combined subsamples with Rinn/r le s exceeding a given 
value yields biased estimators compared to the underlying M-a relation.. The reason for the bias 
in slope is that cuts in R ln fi tend to fall along lines of Mbh oc a 13 ^ 2 (since Ri n a oc Mbhc" 2 and 
Mbh oc a 13 ). This is illustrated by the dashed line of slope 2.0. (Adapted from Giiltckin ct al. 
20093.) 



consists of a sample of 500 galaxies, uniformly distributed in volume out to a distance 
of 30 Mpc. Each galaxy is given a velocity dispersion from a normal distribution in 
log(cr/200 km s _1 ) centered at with standard deviation 0.2. Each galaxy is given a 
BH mass from an M-a relation with intercept a = 8, slope /3 = 4.0, and log-normal 
intrinsic scatter with eo = 0.3 dex. The BH's logarithmic mass is measured with a 
normally distributed measurement error of 0.2 dex and the velocity dispersion has a 
5% error. Since each galaxy has a distance, a BH mass, and a velocity dispersion, we 
calculate -Rinfi = GMbhc -2 and assume the galaxy to be observed with an instrument 
with resolution of d ICS = O'.'l. The results (Fig. [4]) show that cutting out BH masses based 
on their level of resolution would clearly skew fits to the data set since i?i n fl ~ Mbhc -2 
and Mbh ~ a/3 .. so tha t ^infl ~ Mbhc' 3 ~ 2 , resulting in biasing cuts across the data set 



( Gultekin et alj |2009b) 



5. Future Work on Intrinsic Scatter 

The current data set of BH masses is insufficient to test some important questions 
about scaling relations, especially in regards to their intrinsic scatter. One of the most 
important is whether the magnitude o f the intrinsi c scatt er changes across galaxy size. 
This would change the extent of the lLauer et al. ( 2007 ) bias. In order to accurately 
measure the intrinsic scatter across galaxy size, the existing data must be augmented. 
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